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TRACKING OCEAN WAVE SPECTRUM FROM SAR IMAGES 
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ABSTRACT 

For several years we have been working on an end-to-end algorithm 
for recovery of ocean wave spectral peaks from SAR images. We report 
on the progress that has been made, and show that current approaches 
are allowing precisions of 1 percent in wave number, and 0.6 deg in 
direction. 

1. INTRODUCTION 

A SAR impage spectrum is a complex product of the actual ocean 
wave spectrum, the physics of the scattering process, the SAR instru- 
ment, the processing algorithms, and many other effects and param- 
eters. From one or more of these spectra, we seek to either recover 
the actual ocean wave energy spectrum or to estimate other derived 
parameters of interest, such as the locations of storm sources. 

An overall mode- is shown in Fig. 1. The left side shows the 
effects that occur in the production of a SAR spectrum. The right 
side shows the steps we are undertaking in recovery of information 
from such a spectrum. 

Starting from the actual ocean surface, the physics of the micro- 
wave scattering process determines a distribution of scattering cross- 
section on the sea surface, the so-called "oq map". The physics of 
this process was the subject of a workshop held at /'PL in October 1982. 
Wt will not discuss this further here but instead refer the interested 
reader to Beal (1982) . 

The oq map is viewed by the SAP instrument and three types of 
corruption occur: 

1. Dynamic effects: due to motion of the sea surface, scattered 

energy is de focused and/or displaced in the SAR image. Both 
effects serve to decrease azimuth (along-track) resolution, 
creating a roll-off in the system transfer function for 
energy with high azimuth wave number. This is described by 
Monaldo (1983) and Beal, et al.(1983). 

2. SAR instrument: the impulse response of the SAR instrument 

d zermines a modulation transfer function (MTF) , which may 
v&vy with time, and certainly is a function of the SAR pro- 
cessor employed. 
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3. SpecKle: the coherence of the scattering process results in 

constructive or destructive interference of the energy re- 
turned from each resolution cell. This causes the image to 
appear speckled. Viewed as an independent process, the effect 
on image spectra is discussed by Goldfinger (1982) . Below, 
we show what happens when speckling and MTF are considered 
together. 

The final two steps in the production of a SAR spectrum are the 
computations involved in SAR processing and calculation of the image 
spectrum. That SAR processor effects are important is demonstrated by 
Goldfinger (1980) . 

Armed with the "raw" image spectrum, we go through the series of 
steps shown on the right side to produce the final estimate of the 
ocean w **'e energy spectrum or derived parameters. Since the raw spec- 
tra are noisy, it is desirable to average sequential spectra and smooth 
individual spectra. Following this, the system MTF, which has been 
estimated by calibration procedures outlined by Beal, et al. (1981), is 
used to correct for the MTF-speckle effects by an algorithm wa describe 
below. Next, a correction for the dynamic azimuth roll-off is applied, 
as discussed by Monaldo (1983); and further averaging and smoothing are 
performed, if needed. 

Finally, the spectral peaks are located, measured, and tracked 
along the pass. The tracking filter to be used can incorporate inde- 
pendent parameters such as wave height or wind velocity in determina- 
tion of the appropriate filter gains. More will be said of this below. 

If derived data, such as storm source position, are desired, it 
may be necessary to further correct the spectral peak positions for 
systematic effects such as bathymetry or currents before performing 
whatever computations are necessary, such as backward and forward 
extrapolation of *-’ave vectors. 

In presentin' the above model, we have made the cavalier assump- 
tion that the vaiious effects can be separated. For example, we have 
assumed that the dynamic motion effects can be separated from the SAR 
instrument and soeckling effects. Such assumptions should be regarded 
as provisional rather than dogmatic. We expect to test each such 

assumption in turn, and to combine effects when we find they cannot be 

separated . 

For example, in our early work, we assumed that the effects of 
MTF and speckle could be treated independently. However, a theore- 
tical analysis has shown that this i*» not true, and hence the two 
effects are now combined in the model. 

2. PROGRESS ON THE PROBLEM 

For several years, we have been attempting to understand the 

above end-to-end process. Our understanding is still incomplete, but 

partial progress has already been reported in the references. Below, 
we report on the current status of several items. 
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2.1 Scattering Physics 


Although the "transfer function" due to this process is the 
subject of current controversy (Beal, 1982), we have found it useful 
to assume that the relationship between ocean surface slope and radar 
cross-section is linear. We do so for two reasons: 

a. If the transformation between slope and cross-section were 
strongly nonlinear, then the appearance of higher-order 
harmonics and intermodulation products of the fundamental 
ocean frequencies could be expected. Such artifacts have net 
been observed in SAR image spectra gathered in regions of 
open water 

b. The primary nonlinear aspect of SAR ocean imagery does not 
appear to reside in the slope to radar cross-section trans- 
formation but rather in the distortion of the SAR image intro- 
duced by the presence of ocean surface movement. This non- 

li rarity, we have shown, can be modeled as a simple resolu- 
tion loss in a linear system. 

Thus, at present, a linear model is still proving useful for the 
ocean surface slope-radar cross-section relationship. 

2.2 Dynamic Effects 


It is now well established that wave motion at high sea states 
causes an azimuth wave number fall-off that tends to rotate the appar- 
ent direction of the spectral peaks toward the range (cross-track) 
direction (Monaldo, 1983) . Below, we report the progress we have made 
in correcting for this effect. 

2 . 3 SAR Instrument and Speckle 

The effect of coherent speckle on image spectra in a system with 
unity transfer function has been treated (Goldfinger, 1932). If Sq(Jc) 
is the spectrum of the map, the mean SAR image spectrum will be: 


S ' (k) = S Q (k) + 


N h 


VS> 


where N is the i. tber of looks. That is, the mean image spectrum has 
a bias added that is proportional to the integrated spectral power. 

The situation becomes more complex when the system modulation transfer 
function (MTF) is added. The speckle and MTF effects cannot be sepa- 
rated. We have shown that for Gaussian systems: 


M(k) 


s o ( S> 


N J 




where M is the MTF and C is a constant. Thus, the spectrum is now 
weighted by the MTF, as is the power integral taken in the bias term. 
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For non-Gaussian systems, the result is more complex; but analysis of 
the SEASAT MTF has shown that it is well approximated by a Gav^sian, 
so eq. (2) can be used. Analysis cf the statistics of the individual 
spectra shows that the assumptions used in deriving eqs. (1) and (2) 
are valid (Monaldo, 1983) . 

2 . 4 SAR Processing Effects 

In Goldfinger (1982) and (1980) we have investigated the effects 
of processor signature in determining system MTF, and shown how to 
calibrate this aspect of the end-to-end system (Tilley, 1983) . We 
have shown that the total SAR system must be separately calibrated for 
each SAR processor that i-, used. This area is now well enougn under- 
stood that no further woik appears necessary. 

2 . 5 Smoothing and Averagi n g of Spectra 

In the averaging of spectra, we encounter a trade-off between 
averaging more spectra to increase statistical reliability and aver- 
aging fewer spectra to examine changes in the ocean spectrum that 
occur on small spatial scales. Based upon our observation and analy- 
sis of data from pass 1339, which covered a 900 * 100 km region off 
the eastern coast of the United States on September 28, 1978, it wac 
found that spectral variations on the scale of 30 to 40 km tended to 
be random rather than systematic. Evidence based upon measurement of 
the correlation distance of the spectrum tended to confirm that these 
small-scale variations were more characteristic of the instrument and 
scattering processes than of actual variations of the sea surface. 
Accordingly, we chose to smooth spectra with a Gaussian kernal extend- 
ing over about 40 km (i.e., for spectra spaced 6.25 km apart, 15 
spectra were averaged along-track with a Gaussian weighting kernal of 
half width 40 km) . This smoothing appeared to be fairly optimal for 
pass 1339, in that it revealed the local variability caused by the 
Gulf Stream as well as shallow water bottom features. Whether this 
will be generally true for other passes remains to be seen. 

In addition to along-track averaging, we have usually applied 
k-space smoothing to individual spectra. A typical Gaussian smoothing 
distance wan seven wave number intervals. Again a trade-off between 
resolution and statistical reliability occurs. 

2 . 6 Correction for Speckle and MTF 

Based upon eq. (2), we have determined the following algorithm 
for recovering spectra from a SAR image: 

a. Estimate MTF by taking the spectrum of a SAR image of a 
uniform scene (e.g., a field). 

b. Correct the specimen SAR spectra by dividing by MTF and then 
subtracting a bias. 

The efficacy of this algorithm in regions of low to moderate sea state 
has been reported by Beal and Tilley (1980) . The spectral peak 
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positions so obtained have been shown to be highly consistent with 
known storm sources (Monaldo, 1983) . The algorithm is thus well 
established, except that in cases of high sea state it must be com- 
bined with the procedure for azimuth fall-off correction. In general, 
the correction can be applied only out to azimuth wave numbers well 
above the noise. The azimuth fall-off problem apparently can be elim- 
inated only by orbiting very low altitude ('t 250 km) SAR systems. 

2 . 7 Dynamic Corre ction 

We have attempted to correct the azimuth fall-off by using a 
linear semi-empirical model. A Gaussian function in azimuth wave 
number is fit to a region of the spectrum assumed not to contain any 
actual spectral power (i.e., the region of very high wave range number). 
The entire spectrum is then divided by this correction factor. 

When this is done, the position of the observed spectral peaks 
are seen to increase their consistency with known storm sources. Thus, 
empirical evidence exists which confirms the value of this approach. 
However, at this time, the procedure and its limitations are still not 
adequately understood theoretically, a situation we hope to improve in 
the future. 

2.8 Peak Measurement 

So far, three algorithms have been used to determine peak posi- 
tion: 

1. Manual: the spectrum is displayed to an operator who uses 

cross hairs and a joystick to subjectively locate the peak. 

2. Automatic/Peak : the peak value of the spectrum is selected 

automatically. 

3. Automated/Center of Gravity : the center of gravity of all 

pixels exceeding some fraction of the peak is calculated. 

Ideally, when presenting spectra to an operator for manual measurement, 
they should be arranged in random sequence. In our initial study, 
however, this was not done, and the spectra were presented sequential- 
ly. As a result, some "subjective tracking" was performed in obtain- 
ing the manual measurements. 

At the moment, method 3 appears to be easily feasible, and is at 
least as accurate as method 1. Method 2 is most noisy. 

2.9 Tracking 

Most satellite SAR data sets consist of a series of images taken 
at roughly the same time along an extended swath of ocean coverage. 

If the scale of ocean variability is large compared to the separation 
of adjacent images, spectral features can be tracked from one frame to 
the next with improved accuracy of location. 
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To date, we have applied simple smoothing filters and experimented 
with more complex Kalman filters. To construct the latter, we have 
represented the system state vector as the following superposition: 


k(t) = k r (t) + 6k ( t) -> 6k (t) 
- ~u ~s ~c 


(3) 


in which k is the system state, composed of k., the actual wave number, 
6k , small-scale variations, 6k (t) , the occasional jumps due to current 
boundaries. Statistically, we model k 0 a random walk, 6k s as a 
Gaussian process with exponential autocorrelation, and 6k c as a Poisson 
process. 

In choosing filter parameters, we have used the following recently- 
derived information: 

a. On a scale of 10 to 20 km, the variability observed in ocean 
spectra is due to instrumental or noise effects rather than 
systematic variation of the ocean waves. Therefore, tracking 
filters that smooth over such distances will not destroy 
useful information. 

b. The rms errors in spectral peak measurement are dependent upon 
mean wind speed. A simple linear model (c^ * A+BU, where 

is the error in spectral peak measurement an U is wind speed) 
has been fit to the data. Therefore, a filter whose gains are 
dependent on wind speed or other environmental parameters 
would be of value. 

3. RESULTS SO FAR 

Figures 2 and 3 show the results of a typical application of our 
end-to-end estimation procedure to the primary spectral peak of pass 
1339. 


A preliminary analysis of the residuals along the pass shows that 
spectral peaks can be measured with a wave number precision of .0003 
to .0004 rads/m (i.e., 1 percent of mean wave number) and an angular 

precision of 0.6 deg. We emphasize that this is the result of detail- 
ed analysis of a particular wave system over a single pass, and hence, 
should not yet be generalized to other passes or circumstances. 
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